Stable circulation modes in a dual-core matter-wave soliton laser 
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We consider a model of a matter-wave laser generating a periodic array of solitary-wave pulses. 
The system, a general version of which was recently proposed in Ref. 0, is composed of two parallel 
tunnel-coupled cigar-shaped traps (a reservoir and a lasing cavity), solitons being released through 
a valve at one edge of the cavity. We report a stable lasing mode accounted for by circulations of 
a narrow soliton in the cavity, which generates an array of strong pulses (with 10 3 — 10 4 atoms in 
each, the array's duty cycle being ~ 30%) when the soliton periodically hits the valve. 
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*sO ' I. INTRODUCTION AND THE MODEL 

7— I \ 

One of potential applications of Bose-Einstein condensates (BECs) is using them as a source of matter- wave (MW) 
laser beams. Atomic lasing was predicted and demonstrated in a number of theoretical and experimental works, 
starting from early ones PTQ and recently reviewed in Q. Especially interesting is a possibility to design a laser 
operating in a soliton regime, i.e^, generating a sequence of narrow spatial packets of coherent atom waves, as proposed, 
in various settings, in Refs. 0,O|. 

In a recent paper Q , we put forward a model of the soliton MW laser based on a set of two tunnel-coupled parallel 
quasi-one-dimensional (cigar-shaped) traps, one of which is a BEC-filled reservoir, and the other plays the role of a 
I , lasing cavity (the setup is shown in Fig. 1 of Ref. Q), with weakly repulsive and attractive interactions between 
atoms in the reservoir and cavity, respectively (the sign and size of the respective scattering lengths can be controlled 
{SJ [ by external magnetic field through the Feshbach-resonance effect Q ; note that static soliton states in dual-core traps, 
with opposite signs of the scattering lengths in them, were studied in Ref. 0)- In the model, both edges of the 
reservoir, and the left edge of the cavity are impenetrable, while an outlet (valve) at the right edge of the cavity 
releases pulses into an outcoupling MW guide. The model was originally designed with an intention to provide for 
formation of a narrow quasi-soliton that would circulate in the cavity, bouncing from the edges and releasing an 
outgoing pulse each time that it hits the valve. The circulating soliton is supposed to replenish itself by absorbing 
atoms pumped from the reservoir. 

The model is based on a system of linearly coupled Gross-Pitaevskii equations for the wave functions ip(x,t) and 
<p(x,t) of atoms in the cavity and reservoir, with opposite signs in front of the nonlinear terms. The equations are 
cast in a normalized form by scaling h, atomic mass to, and the nonlincarity coefficient in the cavity (in physical 
units, it is 4Trh 2 a/m, with a the negative scattering length of atomic collisions) to be 1, while its counterpart, e, in 
the reservoir is small and positive: 

& : = -(i/2)^xx-HV-^. (i) 

%4h = -(1/2)4^ + E \<f>\ 2 <l>-Kil>, (2) 

with < x < L, i.e., L is the common length of the cavity and reservoir. Keeping e > and the tunnel- coupling 
coefficient n > as free parameters, we use the residual scaling invariance of Eqs. Q and [21 to fix L = 1. The 
model introduced in Ref. [f| also included an external potential and x-dependence of the nonlinearity coefficient 
in equation Q for the cavity; however, the present analysis demonstrates that stable soliton-lasing regimes can be 
achieved without these ingredients. 

As said above, the model assumes that both edges of the reservoir, and the left edge of the cavity are impenetrable, 
the respective boundary conditions (b.c.) being <f>(x — 0) = <f>(x = 1) = ij)(x — 0) = 0. The remaining b.c, 
corresponding to the valve at the cavity's right edge (x = 1) is 

tp x (x = 1) = iqtp(x = 1), (3) 

with q > (as explained in Ref. [jj, it corresponds to the potential drop AC/ = q 2 /2 at the outcoupling edge). After 
straightforward manipulations with Eqs. Q and |2"|). this b.c. yields a balance equation for the total norm of the 
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wave functions ip and <ft, which is proportional to the net number of atoms in the traps, 

N= f (\i>(x)\ 2 + \(j)(x)\ 2 )dx = N 1 +N2. (4) 
Jq 

The balance equation can be cast in the following compact form: 

R(t) = -^-=q\ip(x = l,t)\\ (5) 

Note that the valve gets shut in both limits of q = and q — > oo [in the latter case, Eq. shows that \ijj(x — 1)| 2 ~ 
q~ 2 , hence R vanishes as 1/q]. Monitoring R(t) by means of Eq. (JjjJ, one can characterize the lasing regime (see 
below). The model does not include an explicit form of the outcoupling waveguide (at x > 1), which shapes the 
released pulses into solitons, as soliton formation in a uniform waveguide is a well-studied problem. 

The system is controlled by five parameters: e, k, q, and initial numbers of atoms, and N^°\ in the cavity and 
reservoir, respectively. It is, of course, difficult to exhaustively explore the corresponding parameter space, looking 
for physically meaningful regimes. In particular, a lasing mode based on circulation of a well-defined soliton in the 
cavity was found in Ref. |5j in some parameter region, but it was unstable: after no more than 40 circulations (usually 
fewer), the soliton would quickly come to a halt, getting much broader. Instead of the unstable circulation mode, a 
stable regime of periodic release of pulses was found in a large part of the section of the parameter space explored in 
Ref. [f|. In that regime, a broad "lump" stays put in the cavity, performing cycles of stretching and compression. 
When its right wing periodically reaches the outlet (valve), weak pulses are released. 

The mode of periodic vibrations of the immobile lump is not a spectacular one, and, in practical terms, its drawback 
is that each pulse generated this way contains a small number of atoms, typically not much larger than 10 0. It 
would be more interesting to find a mode of stable soliton circulations in the cavity. In this paper, we demonstrate 
in direct simulations that exploration of a broader parameter region reveals a stable soliton-circulation mode. The 
existence and stability of this mode are directly explained by introducing a simplified analytically tractable version 
of the model, in which we apply the perturbation theory to the soliton circulating in the cavity. The number of 
atoms in each released pulse is A^p U i sc ~ 10 3 — 10 4 , which is quite sufficient to shape the pulses into true MW solitons 
(impulse is larger, roughly, by a factor of 100 than in was provided by the above-mentioned pulsation mode in Ref. 
[5|)- The quality of the soliton array is characterized by its duty cycle (the ratio of the temporal width of each 
pulse to the duration of the cycle), ~ 30%, which makes it possible to avoid the unwanted effect known (in terms 
of telecommunications @) as inter-symbol interference (conspicuous overlap between pulses through their extended 
tails); in the regime investigated in Ref. the duty cycle was ~ 50%. 

The rest of the paper is organized as follows. In Section II, we display generic examples of two different stable 
intracavity-circulation regimes, periodic and quasi-periodic ones, found in direct simulations of Eqs. Q), which 
provide for the MW lasing in the soliton format (we also display an example of another generic regime possible in 
the system, viz., persistent irregular circulations). In Section III, we summarize results of the systematic analysis in 
the form of stability charts for the periodic and quasi-periodic circulation modes in a relevant parameter space, and 
present dependences of essential characteristics of the modes on parameters of the system. The simplified version of 
the model and its analytical perturbative solutions, that demonstrates the existence of the stable-circulation regime, 
are presented in Section IV. The paper is concluded by Section V, where, in particular, we give estimates for physical 
characteristics of the proposed pulsed MW laser. 



II. STABLE CIRCULATION MODES 



Numerical simulations of Eqs. JQ) and J5J were performed by means of a finite-element pseudospectral method. In 
most cases, initial configurations were built (following Ref. (3) as stationary solutions to the same equations, but 
with b.c. replaced by ip(x = 1) = 0. 

Exploring the parameter space of the model in a more systematic way than in Ref. [f| , we were able to find three 
distinct types of regimes featuring persistent circulations of a narrow soliton in the cavity and, accordingly, the MW 
lasing in the pulsed format. The first type is persistent in the sense that the soliton keeps to circulate as long as 
the depletion of the reservoir does not become conspicuous, but the circulations are irregular (quasi-random). An 
example of this regime is shown in Fig. ^ i n terms of the evolution of the density profile in the cavity, \ip(x, i)| 2 , 
and the corresponding release rate, R(t). We do not report more details about this mode, as it seems unusable for 
applications. 

More promising are two other varieties of the circulation mode, one periodic and another one quasi-periodic. 
Typical examples of spatiotemporal patterns specific to these modes are displayed in Figs. |21and|31 by means of both 



3 




(a) (b) 

FIG. 1: A typical example of the spatiotemporal density distribution in the cavity, |^>(a;,i)| 2 (a), and the corresponding release 
rate, defined by Eq. lyj, as a function of time (b), in the case of persistent irregular circulations of the intra-cavity soliton 
circulations. Parameters are k = 0.8, e = 0.0025, q = 32.5, and initial values of the norm of the wave function in the cavity 
and reservoir are 7V( ' = and A^ ' = 1000. 




(a) (b) 

FIG. 2: Generic examples of the spatiotemporal density distribution in the cavity, \ip(x, t)\ 2 , in regimes of periodic (a) and 
quasiperiodic (b) intra-cavity soliton circulations, for = 10, Aj ' = 1000, ft = 0.8, e = 0.002, and (a) q — 12.5, (b) 

q = 21.0. 

spatiotemporal profiles of \ip(x, t)\ 2 and respective contour plots. The corresponding temporal patterns of R(t), which 
determine the profile of the soliton streams released from the cavity (the same way as the temporal shape of an input 
optical signal coupled into a nonlinear fiber determines the form of the soliton array created in it @), are displayed 
in Fig. 01 It is seen that both modes, periodic and quasi-periodic ones, readily establish themselves after a transient 
(this self-establishment was adopted as a criterion of the mode's stability). In the former case, the soliton circulates 
in the central part of the cavity, while in the quasi-periodic mode it sweeps the entire domain and strongly hits the 
valve, in an almost periodic fashion; for this reason, the amplitude of the generated pulses in essentially higher in the 
latter case. Another notable difference of the quasi-periodic regime from its periodic counterpart is period doubling 
obvious in Fig. |3fb) , due to which each cycle of the quasi-periodic oscillations is composed of a narrow tall pulse to 
which a small hump is attached. 

An additional relevant characteristic of the periodic and quasi-periodic regimes is given in terms of the time 
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FIG. 3: The same as in Fig. shown by means of contour plots for the spatiotemporal field of \ip(x, t)\ 2 . 




t t 

(a) (b) 

FIG. 4: The release rate R, defined as per Eq. (SJ, vs. time, in the same cases as in Figs. |5|and|3| Insets in each panel show 
the shape of individual pulses in the established regime, with the continuous lines being a guide to the eye. The duty cycle of 
the (quasi-) periodic pulse arrays, defined as the ratio of the full time width of the temporal pulse at its half-maximum to the 
full period, is 32% (a) and 27% (b). 

dependence of norms N\ and N2 in the cavity and reservoir, as well as the total norm, Ni + N% [see Eq. Q]. For the 
same cases as shown in Figs. [21 and 01 the dependences are displayed in Fig. [5] They clearly highlight the character 
of the dynamical regimes: a relatively fast matter exchange between the reservoir and cavity, and slower depletion of 
the total norm, which actually happens by jumps when a new pulse is released into the outcoupling MW waveguide. 

All stable periodic and quasi-periodic soliton-lasing modes found at other values of parameters (see stability charts 
presented in the next section) are very similar to the examples displayed in Figs. [5] -03 As concerns the persistent 
irregular regimes, they are quite similar to the example displayed in Fig. 2J. 
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(a) 



(b) 



FIG. 5: The evolution of the total norm, Ni + N2 (rhombuses), and norms in the lasing cavity and reservoir, iVi (squares) and 
N2 (triangles), in the same established dynamical regimes [periodic and quasiperiodic ones, (a) and (b)] which are shown in 
Figs. Hand 01 

III. STABILITY MAPS FOR THE CIRCULATION MODES. 

As the model contains several parameters, it is necessary to identify those which strongly affect the stability of the 
lasing modes. For this purpose, we summarize results of a large number of simulations in charts displayed in Figs. 
El and The periodic and quasi-periodic circulation modes are stable, respectively, inside the area bounded by the 
continuous line, and between dashed and continuous borders. The charts displayed in these figures were obtained by 
varying q and e in steps of 2 and 0.001, respectively, unless strong sensitivity to variation of the parameters dictated 
to use reduced steps. 

Comparison of Figs. OJa) and|^b), which differ by values of the initial number of atoms in the cavity, demonstrates 
that the outer stability border (appertaining to the quasi-periodic regime) is relatively insensitive to this parameter, 
while the inner border, which separates the stability areas of the periodic and quasi-periodic regimes, is conspicuously 
affected by the initial conditions. The effect of variation of the strength of the linear coupling between the reservoir 
and cavity, k, on the stability region is illustrated by comparison of Figs. Elb) and [7] Naturally, the decrease of k 
entails shrinkage of the stability area, as it becomes harder for the circulating soliton to replenish itself by collecting 
matter pumped from the reservoir (see also the analytical model presented in Section V). In addition, the results 
suggest that, at smaller k, the stability region shifts to higher values of q. Detailed analysis shows that the increase 
of q makes the release rate R smaller in this case, thus compensating the decrease of the pumping rate due to the 
reduction in the coupling strength. 

A noteworthy peculiarity of Fig. is a nearly flat lower stability border, which is found at q = q m i n ~ 4. An 
explanation for this flat border is provided by the analytical model presented below. 

Beneath lower borders of the stability areas in Figs. 4 and 5, no circulation regime is possible. Above the upper 
stability borders, persistent irregular circulations take place, similar to the example shown in Fig. |3 Outside the 
stability area, there is an additional border between the no-circulation and irregular-circulation regimes, which we do 
not display here, as it is of little relevance to applications. 

Effects of the variation of control parameters on characteristics of the generated pulse arrays were investigated too. 
In particular, a change in does not affect the results in any tangible way, and arrays of the outcoupling pulses are 
quite insensitive to the variation of e too. On the other hand, the dependence on the valve parameter q is conspicuous, 
as seen in Fig. 



IV. THE ANALYTICAL MODEL 



A qualitative explanation to the stable regime of periodic circulation of the soliton in the cavity may be provided 
by the perturbation theory, which, at the zero order, assumes shuttle oscillations of a soliton in the cavity of length 
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(b) 



FIG. 6: Stability borders of the periodic and quasi-periodic regimes in the parameter plane (e,q), for k = 0.8, iVj = 1000 
and (a) N[ 0) = 0, (b) N[ 0) = 10. 
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sility border of quasi-periodic mode 
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FIG. 7: The same as in Fig. 4(b), except that k = 0.5. 



L with zero b.c, ip(x — 0) = ip(x = L) = 0, with no coupling to the reservoir, n = 0. The respective zero-order 
approximation for the soliton is 



ip Q (x,t) = »/sech(»7(z-£(t)) e mx+l * {x \ 
d£, d$ 1 / 2 a\ 



(6) 



where 77 and v are its amplitude and instantaneous velocity. In this approximation, Eq. (JQ conserves two dynamical 
invariants, viz., the norm, N% = J Q \ip(x)\ 2 dx, see Eq. Q (in this section, we do not set L = 1, to present the 
dependence of the results on parameters in a clearer form), and the energy (Hamiltonian), 



^1 = 7^ 



(|^| 2 -M 4 ) dx. 



(7) 



Note that, for unperturbed soliton (jfJJl, the values of these invariants are 



(Ni) sol = 277, (iJr) sol = V v 2 - -n 



(8) 



Zero b.c. imply the presence of formal mirror images of the soliton beyond the points of x = and x — L, i.e., 
il>(— x) = —ip{x), and tp(x) — —ijj(2L — x), respectively. As is well known, a pair of solitons with opposite signs repel 
each other |!j , which explains bounces of the soliton from edges of the cavity and its shuttle motion. 
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FIG. 8: The maximum value ("amplitude") of the release rate R, and its oscillation period, versus the outcoupling coefficient 
q, for k — 0.8, e = 0.002, = 1000, and two different values of n[°\ Data sets for n[ ^ = terminate at points beyond 

which stable circulation modes are not found. 



In the first nontrivial approximation, we treat the coupling coefficient k and b.c. coefficient q in Eq. JSJ) as small 
parameters. To predict the possibility of stable periodic circulations of the soliton in this case, we will treat it as a 
quasiparticle trapped in an infinitely deep potential box, between x = and x — L. Then, the shuttle motion can be 
predicted by means of balance equations for the norm and energy [this approach was elaborated in detail for a kink 
(topological solitons) moving as a shuttle in models based on a perturbed sine-Gordon equation 

First, note that Eq. J^J makes it possible to calculate the loss of norm Ni for the soliton bouncing from the right 
edge (x = L): 



/+oo 
\^{x = L,t)\ 2 dt = -2^-, 
-OO ^ 



(9) 



where ip{x = L) was taken as per the unperturbed soliton 0. As concerns the Hamiltonian, in the absence of 
perturbations its conservation follows from the continuity equation relating the Hamiltonian density and current 
conjugate to it, 



d_ 

dt 



l# 



Integration of this equation over the domain of < x < L gives rise to the evolution equation for Hi similar to its 
counterpart © for the norm, 



dHi 
~dt 



= q 



(10) 



Substituting unperturbed soliton © in this relation and subsequent integration yield the net change of Hi induced 
by the bounce of the soliton from the valve, cf . Eq. (JHJ : 



(A#i) c 



+ 00 



dHi qrj , 2 2 ^ 
dt v ' 



(11) 



The change in the norm and energy incurred by the collision of the soliton with the valve is to be compensated by 
matter which was pumped from the reservoir and absorbed by the soliton sweeping the cavity. Assuming, in the spirit 
of the perturbation theory, that the matter-wave background built up in the cavity between passages of the soliton, 
ipo(x), has a small amplitude, and the soliton collects all the available matter, the linearization of Eq. JIJ yields 



ip Q (x) = incj)o(x)T(x), 

where r(x) is the time between two passages of the soliton through a given point, and 



[2N^ . (irx\ 

\j — sm {T) 



(12) 



(13) 
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is the ground-state wave function in the reservoir, with the number of atoms N 2 in it; the depletion of the reservoir 
is neglected, which is a natural assumption in the perturbative setting. As the soliton passes each point twice in the 
course of the cycle of circulations, it is easy to see that t(x) has two branches, 

ti = 2x/v, t 2 = 2 (L - x) /v. (14) 

Then, it follows from Eqs. (@J, J7J and (|14|) that the norm and energy collected by the soliton within the circulation 
cycle are 

,2 r 2 / \ , 2 / \i 7 ^ k 2 N 2 L 2 



(ANi) K = (Kfrr I [rl{x) + r 2 {x)] dx = C N —^—, (15) 

(k 2 N 2 ) 2 L 3 



(AH 1 ) K = -i (^o) 4 ^ [rt(x) + ri(x)] dx = -C H - 



(16) 



where C N = 4 (| - « 0.93 and C H = f - f + « 2.25. 

The periodic circulations are provided by obvious balance conditions, (ANx) q + (AN X ) K = and (AH X ) +(AHi) K — 
0. Using Eqs. Ijlljl and (|15|l . 1(16() . these conditions determine equilibrium values of the soliton's amplitude and 
velocity, 770 and vq- 



, C N k 2 N 2 L 



= -f—^-, (17) 
_C H (^iV 2 ) 2 L3 

wo = Vo -0v , ov = — c ■ (18) 

2 577° 

Within the framework of the analytical model it is also possible to analyze stability of this equilibrium solution, by 
assuming small deviations of r\ and v from values (|17|) and l|18|l . For the corresponding nonequilibrium solutions, Eqs. 
©, CU and O, (EH) give rise to a map, 2?? -> 2rf = 2r]+(AN 1 ) q + (AN 1 ) K , H X ^H[ = H x + (AJ?i) g + {AH X ) K [in 
fact, it is more convenient to analyze the map written in terms of ij and reduced Hamiltonian, h = Hi — (2/3)r/ 3 ]. Then, 
equilibrium solution $Tj$ and 1)18(1 is a fixed point (FP) of the map, and its stability is determined by multiplicators 
/i, i.e., eigenvalues of the map linearized around the FP. Direct calculation yields two eigenvalues, (i v = 1 — 2 (q/r)o) 
and /i/j = 1 — {q/rjo). As g is a small positive parameter, we conclude that < Hn,^h < 1, hence the FP is stable 
indeed. 

Despite many approximations adopted in the analytical version of the model, the prediction for the soliton's 
amplitude given by Eq. (|17f) generally agrees with numerical results. For instance, in the example of the periodic regime 

displayed in Figs. IHa),|3{a) and^Ja), the averaged maximum density of the solitonic waveform is / ( \ip(x, t\ 2 ) \ ~ 

\ V / max/ 

27, while Eq. fTTjl predicts, for the same case, r] 2 « 25. 

The analytical model also helps to explain the presence of the lower border of the stability region in Fig. HJ1 Indeed, 
stable circulations of the soliton are possible if the release of matter through the outlet can compensate the permanent 
inflow of matter from the reservoir. Adopting a rough estimate according to which 77, v and q are on the same order 
of magnitude, which is borne out by numerical data in the general case, and making use of Eqs. © and JTHJ, one 
can put the latter condition in the form of g 3 > q^ in ~ k 2 N 2 L 2 . For parameter values corresponding to Fig. this 
yields g m ; n ~ 7, which generally agrees with the situation shown in that figure. 

One can also understand why the lower stability border in Fig. is nearly flat (in fact, the explanation is quite 
a general one, and is not predicated on the validity of the above analytical model) . Indeed, it follows from Eq. (J2J 
that nonlinear deformation of ground state 1|13[) is negligible under condition e <C e m ax = 71-2 / (N 2 L). In the case 
corresponding to Fig. |B| with N 2 = 1000 and L = 1, one obtains £ max — 0.01, which complies with what one can 
observe in Fig. 

A limitation of the quantitative validity of results produced by the analytical model can be seen from the calculation 
of the maximum value of the matter-release rate: Eqs. (JSJ, 10 and l|17fl yield a perturbation-theory result for it, 
-Rmax = (Cjv/2) k 2 N 2 L 2 , which does not depend on the outlet parameter, q. On the other hand, Fig. [S]shows that, 
in the range of 5 < q < 17.5, the numerical results yield an approximate dependence, i? max ~ g~ ' 75 . However, the 
perturbation theory does not really apply in this range. 



V. CONCLUSION 



In this paper, we have investigated the model of the matter- wave laser generating a periodic array of coherent pulses, 
which was recently introduced in Ref . . While only a relatively primitive regime of vibrations of a broad immobile 
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lump in the cavity was found in the original consideration of the model, here we have demonstrated that exploration 
of a broader parameter region reveals stable lasing modes accounted for by periodic or quasi-periodic circulations of 
a narrow stable soliton in the cavity; the soliton periodically generates an outcoupling pulse, hitting the valve, and 
then gets replenished, absorbing matter pumped from the reservoir. The existence of the stable circulation regime, 
and some features of the numerically found stability regions were explained in the framework of a simplified analytical 
model, that treated the intra-cavity soliton as a quasi-particle, and predicted the stable regime by means of balance 
conditions for the number of atoms and energy in the cavity. 

Getting back from the Gross-Pitaevskii equations in the normalized form of Eqs. J3J) to physical units, and 
making use of typical values of physical parameters for the 7 Li condensate, it is easy to conclude that the circulation 
regime found in this work easily provides for the generation of regular pulses built of J\f = 10 3 — 10 4 atoms each (the 
vibrational lasing mode investigated in Ref. p| gave rise to Af smaller by a factor of ~ 100). With the scaled value of 

the initial norm, = 1000, employed in the examples, the laser would generate ~ 100 pulses before the reservoir 
depletion will becomes tangible. However, the latter number may be readily made indefinitely large, as indefinite 

increase of does not affect the stability of the circulation regime. 

Undoing the rescaling which cast the Gross-Pitaevskii equation in the normalized form of Eq. Q , it is easy to see 
that, for an experimentally relevant value of the cavity length, L — 100 /im, the time unit in the normalized equation 
corresponds to physical time to = mL 2 /h w 1 s (recall m is the atomic mass, i.e., m w 1.2 x 10~ 26 kg for 7 Li). 
Therefore, according to Fig. 03 a typical circulation period is ~ 0.2 s in physical units. Lastly, detailed consideration 
of numerical data shows that the typical value of the duty cycle characterizing the quality of the generated pulse 
array is ~ 30% (the quality of small-amplitude arrays found in Ref. Q was poorer, with the duty cycle ~ 50%, i.e., 
the pulses were more overlapped than in the circulation regime). 

The realization of the proposed matter-wave soliton laser by means of available experimental techniques seems to 
be quite feasible. It may also be possible to consider another variant of the proposed matter-wave soliton laser, with 
outlets at both edges at the cavity. The so modified system may generate a double stream of matter-wave solitons, 
which, however, needs separate consideration. 

Acknowledgement 

The work of B.A.M. was supported, in a part, by the Israel Science Foundation through the Center-of-Excellence 
grant No. 8006/03. 



10 



[1] Holland M, Burnett K, Gardiner C, Cirac J I, and Zoller P 1996 Phys. Rev. A 54 R1757; Moy G M and Savage C M 1997 
ibid. 56 R1087; Kneer B, Wong T, Vogel K, Schleich W P, and Walls D F 1998 ibid 58 4841; Band Y B, Julienne P S and 
Trippenbach M 1999 ibid. 59 3823; Martin J L, McKenzie C R, Thomas N R, Sharpe J C, Warrington D M, Manson P J, 
Sandle W J and Wilson A C 1999 J. Phys. B - At. Mol. Opt. Phys. 32 3065 

[2] Mewes M O, Andrews M R, Kurn D M, Durfee D S, Townsend C G, and Ketterle W 1997 Phys. Rev. Lett. 78 582 Miesner 
H J, Stamper-Kurn D M, Andrews M R, Durfee D S, Inouye S and Ketterle W 1998 Science 279 1005; Bloch I, Hansen T 
W and Esslinger T 1999 Phys. Rev. Lett. 82 3008; Deng L, Hagley E W, Wen J, Trippenbach M, Band Y, Julienne P S, 
Simsarian J E, Helmerson K, Rolston S L and Phillips W D 1999 Nature 398 218; Hagley E W, Deng L, Kozuma M, Wen 
J, Helmerson K, Rolston S L and Phillips W D 1999 Science 283 1706 

[3] Bongs K and Sengstock K 2004 Rep. Progr. Phys. 67 907 

[4] Leboeuf P, Pavloff N, and Sinha S 2003 Phys. Rev. A 68 063608; Carr L D and Brand J 2004 Phys. Rev. A 70 033607; 

Rodas- Verde M I, Michinel H and Perez-Garci'a V M 2005 Phys. Rev. Lett. 95 153903 
[5] Chen P Y P and Malomed B A 2005 J. Phys. B - At. Mol. Opt. Phys. 38 4221 
[6] Agrawal G P, Fiber-Optic Communication Systems (John Wiley & Sons, Inc.: New York, 1997) 

[7] Kagan Y , Surkov E L, and Shlyapnikov G V 1997 Phys. Rev. Lett. 79 2604; Roberts J L, Claussen N R, Burke J P, 
Jr, Greene C H, Cornell E A, and Wieman C E 1998 ibid. 81 5109; Inouye S, Andrews M R, Stenger J, Miesner H-J, 
Stamper-Kurn D M, and Ketterle W 1998 Nature 392 151 

[8] Shchesnovich V S, Malomed B A, and Kraenkel R A 2004 Physica D 188 213; Shchesnovich V S, Cavalcanti S B, and 
Kraenkel R A 2004 Phys. Rev. A 69 033609 

[9] Kivshar Y S and Malomed B A 1989 Rev. Mod. Phys. 61 763 



